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Abstract 

p^i . We show that the concept of complete symmetry group introduced by Krause (J. Math. 

>- \ Phys. 35 (1994), 5734-5748) in the context of the Newtonian Kepler problem has 

wider applicability, extending to the relativistic context of the Einstein equations de- 
scribing spherically symmetric bodies with certain conformal Killing symmetries. We 
also provide a simple demonstration of the nonuniqueness of the complete symmetry 

O '. group. 

o. 

*£ . 1 Introduction 



In 1994 Krause [12| introduced the concept of a complete symmetry group in classical me- 
J> \ chanics. This concept was somewhat stronger than that of the standard symmetry group 

generally used at the time, which was the group represented by the Lie point symme- 
tries of the system of differential equations describing the dynamical system. A complete 
symmetry group realisation in mechanics was required to be endowed with the properties 
that (i) the group act freely and transitively on the manifold of all allowed motions of the 
system and (ii) the given equations of motion be the only ordinary differential equations 
that remain invariant under the specified action of the group. The vehicle which Krause 
used to demonstrate the concept was the classical Kepler problem described in reduced 
coordinates by the second order ordinary differential equation 

r + ^ = 0, (1.1) 
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where \x is a constant, r = (xi, x 2 ,x 3 ) and r = \r\ in a standard notation. The Lie point 
symmetries of the system are given by 

G\ = x 2 d X3 - x 3 d X2 , G 2 = x 3 d Xl - x x d X3 , G 3 = x\d X2 - x 2 d Zl , 

G A = d t , G 5 = 3td t + 2rd r , (1.2) 

where G%, G 2 and G 3 are the generators of rotation, G4 is the generator of time translations 
and G5 represents rescaling which maps orbits of a given eccentricity into orbits of the same 
eccentricity and is the basis of Kepler's Third Law because of the happy happenstance 
that for the Kepler problem (and also the simple harmonic oscillator) the eccentricity 
does not enter into the relationship between the period of the orbit and the length of the 
semimajor axis as it does for other periodic orbits || . The algebra of the point symmetries 
is so(3) © A 2 . In addition to the standard Lie point symmetries Krause determined the 
three nonlocal symmetries 



X\ = ( 2 / x±dt j d t + xird r , X 2 = I 2 / x 2 dt j d t + x 2 rd r , 

X 3 = ( 2 / x 3 dt J d t + x 3 rd r (1.3) 



and with these eight symmetries was able to specify the precise form of ( |1 . 1| ) up to the 
value of the constant \i. That the constant fi could not be specified is not surprising since 



it can be removed from the equation (1.1) by a rescaling of the vector r. 

The work of Krause was and is very interesting. Apart from any other features which 
are peculiar to the Kepler problem, there was presented the new approach to the calcula- 
tion of nonlocal symmetries which are not easy to calculate in an algorithmic way m . It 
is unfortunate that in his paper Krause made the comment that the complete symmetry 
group could not be obtained for the Kepler problem by the standard analysis used for 



Lie point symmetries. Rising to the implied challenge Nucci [23] proceeded to devise a 
scheme of reduction of order of the system, in essence a removal of time as the indepen- 
dent variable so that the coefficient function of dt appeared only as the derivative and 
hence a local function for the type of nonlocal symmetry used by Kause, which made 
the determination of the nonlocal symmetries above a process of the determination of the 
Lie point symmetries of the reduced system. More recently Nucci and Leach |24J have 
shown that the number of nonlocal symmetries of this class for the Kepler problem was 
somewhat greater than Krause had reported. Furthermore, in contrast to his statement 
that one does not obtain a Lie algebra when the nonlocal symmetries are considered, they 
presented the algebra. 

Krause emphasised that his treatment applied to mechanical systems based on the New- 
tonian laws, but there can be no doubt that the concept of a complete symmetry group is 
applicable to all systems of ordinary differential equations which possess a suitable number 
of symmetries to define precisely the system of equations. Krause further remarks that it 
is not known whether the specific realisation of a complete symmetry group exists for any 
given Newtonian system -- more generally system of ordinary differential equations - 
and that the question was open. There is a specific realisation for one-dimensional linear 
second order equations, sl(3, R), but even for nonlinear one-dimensional systems the prob- 
lem was unresolved. This immediately raises the question of what the situation would be 
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in the case of dynamical systems known to display the characteristics of nonintegrability. 
Considering that the existence of symmetry is regarded as one of the criteria relating to 
integrability, the concept of the complete symmetry group of a nonintegrable dynamical 
system would appear to present something of a paradox. 



In this paper, which is the fifth of a series |16|, [L7|, [Ll], |2^] devoted to investigating the 
connections between the three main topics in dynamics, videlicet symmetry, singularities 
and integrability, topics which superficially are unrelated in their mathematics and yet are 
intimately intertwined, we demonstrate that it is possible to have a complete symmetry 
group of an ordinary differential equation which does not have sufficient Lie point sym- 
metries for integrability and is not integrable in the sense of Painleve. In the process we 
demonstrate the resolution of the apparent paradox. 

In the normal course of events one associates the existence of symmetry in a differen- 
tial equation (or system of differential equations) with the integrability of the equation. 
However, the normal course of events also restricts the consideration of a symmetry to the 
class of point, contact or generalised symmetries. In the very introduction of the concept 
of a complete symmetry group Krause was obliged to resort to nonlocal symmetries for 
the complete specification of the system of differential equations for the Kepler problem 
which is well-known to be an integrable system. The need to use nonlocal symmetries 
for the complete symmetry group is not automatic. For example the equation y" = is 
specified completely by three point symmetries. 

In this paper we are not concerned with the existence of symmetries which will enable 
the equation to be integrated. We are concerned with the existence of symmetries which 
will enable the equation to be specified completely, perhaps up to an arbitrary parameter 
as was the case of Krause's complete symmetry group for the Kepler problem. The role 
played by the elements of the complete symmetry group is to fix the form of the dynam- 
ical equations. Our aim is to demonstrate that this be the case whether the system be 
integrable or not integrable. As we shall see, the very nature of the symmetries obtained is 
not likely to cause one to imput integrability of the model equation. It would appear that 
symmetries can play two roles, one to specify the nature of the equation and the other to 
provide part of the route to integrability. A priori one cannot expect both roles to the 
played by the one symmetry. 

This paper is structured as follows. In the next section we review the known properties 
of our model equation. In Section 3 we calculate its symmetries and demonstrate that 
within the number of these symmetries there is a subset which is the complete symmetry 
group of the model equation. In Section 4 we conclude with some observations. 

2 Painleve nonintegrability 

Our model equation is 

y'" + y" + yy' = o, (2.1) 

where the prime indicates differentiation with respect to the independent variable x. Equa- 
tion ( |2.ip arises in general relativity. On the assumption of a conformal Killing vector of 
a particular type for a spherically symmetric shear-free gravitational field Dyer, McVittie 
and Oattes Q obtained the third order field equation 

V 2 T^ + fi(2m - 1)T W + (m 2 - 2m + 2T) T M = 0. (2.2) 
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Under the transformation 

/i = exp f — ^— r j , T(^) = 2(m-2) 2 y(x)-Um 2 -Am + 3) (2.3) 

V 2m — 4 / 2 



we obtain (2.1) from (2^). To determine the leading order behaviour of (2A) we make the 

substitution y = a\ p ', where x = x ~ x o into ( |2.1| ), We find that the first and the third 
terms are dominant, a = — 12 and p = — 2. The substitution of 

y = -12 X " 2 + p X r ~ 2 (2.4) 

into the two dominant terms gives the resonances r = —1,4, 6. We substitute 

4 

y = -12y~ 2 + Y. a ^ ( 2 - 5 ) 

into the full equation ( |2.lD to see it is compatible at the resonances. We find that 

12 1 1 . N 

a - 1 = Y' ao= 25 &nd ai= 125 (2 - 6) 

for the coefficients up to the first resonance at r = 4. At the first resonance we have 
compatibility and, by construction, a<i is arbitrary. The next coefficient is given by 

1 + 1250a 2 ,__. 

a-i = . (2.7) 

S 18750 V ; 

However, at the resonance, r = 6, the condition of compatibility requires that 

1 , s 

a 2 = , (2.8) 

1250' v ; 



i.e. the coefficient 03 is zero. Consequently (2.1) fails the Painleve test. The solution 



requires logarithmic terms to be introduced into the Laurent expansion and this means that 
the general solution has infinitely many branches. Infinite branching is strongly connected 
with nonintegrability and this suggests that (2.1) could exhibit chaos [|l|][|28|, p. 348][]. 



Analysis of the numerical properties of the equation reveals that there is instability in the 
numerical solution but not chaos. 

Even a cursory examination of ( |2.lD reveals the existence of a first integral 

i(y,y\y") = y" + y' + \y 2 (2.9) 

which is associated with the obvious Lie point symmetry of ( |2.1[ ), videlicet d x . Instead of 
viewing (|2.9|) as a first integral of fl2.1|) one can consider it as an equation of the second 
order containing a parameter, /, videlicet 

y" + y' + \y 2 -I = U. (2.10) 



lr rhe logarithmic singularity of the general solution can be avoided by accepting the constraint, (|2 

5 A T-Sd-r 125 A + 1250 



on the coefficient ai so that one obtains the series expansion, y — —12\ 2 + ^ X X + 25 + t|fX + T^TTX^ + 



<24X 4 + • ■ ■ • This expansion contains only the two free parameters, xo and 04 and is interpreted to be the 
local representation of a particular solution of (EJ) in the sense of Cotsakis and Leach 3, uM. There is 
another interpretation which is that this series expansion is the ge nera l solution of a second order equation 
which, in some sense, is an invariant of the third order equation, (|2.r|). 
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This equation can be written as a family of two-dimensional first order systems, 

x = y, (2.11) 

y = -y- l -x 2 + I. (2.12) 

The equilibrium points of this system occur only for I > 0. They are x = ±V2I, y = 0. 
A simple application of the Bendixon-Dulac criterion reveals that there are no limit cycles. 
Thus, as I goes from negative to positive values, the system experiences a supercritical 
bifurcation with an additional qualitative change in the nature of the critical points. The 
equilibrium points on the positive x axis are always saddle points. On the other hand 
those lying on the negative x axis are stable nodes for / < 1/32 and stable foci with 
clockwise direction for I > 1/32. 

It is sometimes natural to consider the dynamics of a two-dimensional family of systems 
depending upon parameters as equivalent to that of a higher order model equation from 
which it can derived. In our case, however, since the two-dimensional reduced system has 
bifurcations and lines of equilibria it may have completely different (and simpler) dynamics 
than that hidden in the third order equation. In fact we can apply the Painleve test to 
( 2.1 0| ) and find that it satisfies the necessary conditions for the possession of the Painleve 



property only in the case that / = —18/625. Then it has the expansion 

12 111 

V = -12y" 2 + — Y + 25 + Y X 2 X 3 + hx 4 + ■■■ , (2-13) 

y a -r 5 a 125 a 1250 o A 187500 A A y ' 

where 64 is the arbitrary constant introduced at the resonance r = 6. The integrability of 

y" + y' + I y 2 + 2i =0 ( 214) 

y y 2 y 625 v ' 



has been demonstrated by Ince [1C] who transformed ( 2.14| ) to 

d 2 w 



dz 2 W, (2.15) 

the solution of which can be expressed in terms of elliptic functions, by means of the 
Kummer-Liouville transformation |13|, 20] 

/ 2z\ 6 hi f z\ 

y(x) = expi- — Jw{x) + —, x = — -=exp\^--j. (2.16) 



Equation (p. 15) is integrable in the sense of Painleve and the transformation (2.16J) is 



single- valued. Hence ( 2.14| ) is integrable. Whatever happens to the solutions of (|2.1| 



which commence at some point not on the surface defined by ( |2,10D in the phase space 
with the particular value of / = —18/625 is not obvious, but those which do are regular 
and remain on the surface. 

As a final note on the properties of our model equation we remark that the equation, 
(|2.14| ), has the additional Lie point symmetry (calculated using Program Lie (^, |27|), 
exp (I) \d x — (|y — y^) d y ] , in addition to the symmetry, d x , possessed for general values 
of the parameter, /. The algebra of the symmetries is Type III of Lie's classification of 
two-dimensional algebras [19, p. 430, Kap. 19] which has the normal form 



G 1 = d y , G 2 = xd x + yd y (2.17) 
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so that [Gi,Gr2] = G\ and the normal form of the equation invariant under the two 
symmetries ( [2.17 ) is 

XY" = F(Y'). (2.18) 

The tidy form of the equation, ( |2.15| ), arises from using the nonnormal form 

Gi = d x , G 2 = -xd x + yd y . (2.19) 



The second symmetry in (2.17) should not be expected to be related to the complete 
symmetry group to be computed since ( |2.14| ) is not equivalent to (|2.1| ) as the particular 
value of / brings us into the realm of configurational invariants ||, |2(| . 

3 Lie symmetries 



Equation ( |2.lD has just the one Lie point symmetry, d x , which cannot be the complete 
symmetry group of (^4]) since the equation 

y'" = f{y,y',y") (3.1) 



is the general form of a third order ordinary differential equation invariant under 
Examination of (^4]) for contact symmetries shows that there are none. Consequently the 
complete symmetry group (we now speak strictly in the sense of Krause) must comprise 
nonlocal or generalised symmetries. To calculate these we follow the procedure outlined 
in Pillay et dl [pH] and used with great effect by Bouquet et al |l5[ in their study of second 
order ordinary differential equations possessing integrating factors j|] . 

When one studies symmetries of differential equations without placing any restrictions 
upon the coefficient functions one may replace the differential operator 

G = i8 x + rjd y (3.2) 

with the equivalent operator 

G = fjd y , (3.3) 



where rj = r\ — y'£. The advantage of using ( |3.3| ) is that the calculations are simpler. 
Dropping the overbar we find that the action of the third extension of (|3.3[) on (£4]) gives 



rf" + v" + v'y + w' = 0- (3.4) 

Our task is to find a sufficient number of independent functions, 77, which are solutions 



of (|3.4|) to be able to define (|2.l|) uniquely and so constitute the required number of 
symmetries to make a complete symmetry group. 

It is a trivial matter to integrate ( [3,4| ) once. We have 

r]" + v ' + r)y = C, (3.5) 



2 One could consider introducing the idea of degrees of completeness in symmetry groups. This would 
depart from the criterion of Krause that the equation be unique. However, there is no reason why we 
should not elaborate on Krause's proposition. 
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where C is an arbitrary constant of integration. We replace the coefficient of r/ from the 
original differential equation, ( |2.lD , to obtain 

(3.6) 



(v" + v')y'-v(y'" + y") = c y ' 

which can be integrated by parts to give 



rjy' - rjy" + / (r/'y' - r]y")dx = Cy + K, 



(3.7) 



where K is another arbitrary constant. If we let uj = J(rj'y' — rjy")dx, the integro- 
differential equation (|3.7|) becomes the first order linear nonhomogeneous equation 



J + & = Cy + K 



which has the solution 



uj = Je x + e 



X I _— X 



J{Cy + K) 



e x dx, 



(3.8) 



(3.9) 



where J is a constant of integration, so that r\ now satisfies its own linear first order 
nonhomogeneous equation, 



v 'y' - V y" = Cy + K - Je~ x - e~ x I ' (Cy + K)e x dx. 
The solution of ( 3.1C| ) is 

e~ x f(Cy + K)e x 



ri = Ay' + y' j -^ 



Cy + K - Je~ 



dx 



dx, 



(3.10) 



(3.11) 



where A is a constant of integration. 

In ( 3.1 1| ) there are four constants, but one of them, K, is spurious. We have the three 
symmetries 

G 1 = y'd y «• Gi = 8 X , 

G3 = y'(J^(y-e- x Jye x dx))d y 

** G z=(Jy^(y- e ~ X Jye X ^d x , (3.12) 

where we have added the corresponding operator when we use only d x for the symmetry. 
We observe that the first of ( 3,12; ) is just the obvious symmetry of (2.1). 

Our immediate task is to find the form of the third order differential equation which 
is invariant under the actions of the third extensions of these three symmetries. We 
commence with the general equation 



y'" = f(x,y,y',y"). 



(3.13) 
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For the calculations we use the d x form of the symmetries. Invariance under G\ implies 
that the equation have the form 

y'" = f(y,y',y"). (3.14) 

Invariance under G2 and G3 is not obtainable by inspection! The third extension of G% is 



G. 



[3] 



y,2 eX l 9 - 



d v > + I -r-r ) <9„" 



y > e x 1 v 



y'e x I v 



y'" + y" + y' 

y'2 e X 



On 



(3.15) 



and, when it acts on ( 3.14 ), we obtain the linear first order partial differential equation 
which has the two characteristics 



u = y" + y and 



/ + y' + y" 
y' 



so that the right side of ( |3.14| ) has the form 

f = -(y" + y') + y'g{y" + y',y), 

where g is an arbitrary function of its arguments. 
The third extension of G3 is 



(3.17) 



(3.18) 



G 



[3] 



dx 

7/2 



y-e 



/ ye x dx ))d x - — (y-e x ye x dx J d y 



7 I y -y + e" 



/ ye x dx ) dy// + d y m . 



(3.19) 



From the action of ( p.!9| ) and (3.14) with / as given in ( |3.18| ) we find that g is inde- 
pendent of its first argument and so the general form of the third order equation invariant 
under the three symmetries (3.12) is 

y"' + y" + y' + y'Ky) = o, (3.20) 

where h(y) is an arbitrary function. 

It is quite evident from ( |3.20| ) that the three symmetries ( |3.12| ), which arise from 
the direct integration of the linear third order ordinary differential equation (|3.4|), are 



insufficient to specify completely our model equation, (24), and so we must find some 
other symmetry. We assume that the symmetry has the form S,d x . The action of the third 
extension of the symmetry on ([2.1D gives 



y'e + (V + y') f + (2y'" + y") t = o, (3.21) 

in which (|2.1| ) has been used. Equation (|3.21[ ) is exact and can be integrated once to give 

y>Z» + (2y" + y ')e = A, (3.22) 
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where A is a constant of integration, which is a linear first order nonhomogeneous equation 
for £' with the solution 



?=a 



f y'e x dx 1 



y' 2 e x 



+ B 



y' 2 e x 



and so 



£ = A 



J y'e x dx 
y' 2 e x 



dx + B 



dx | c 
y ;2 e x 



(3.23) 



(3.24) 



where B and C are further constants of integration. Unfortunately, and this should not 
be surprising, the three symmetries in (3.24) are just the three symmetries listed in ( p. 12 ). 
Evidently we are not going to be able to obtain the necessary additional symmetries 
the easy way and so we must assume that the symmetry has the form (c/ the Ansatz of 
Krause @) 

G=(hdx\d x + r,d y . (3.25) 

The action of the third extension of ( p. 25 ) on the equation (|2.l] ) gives 

y'i" + y"i' + (V + y')C + ( V" + y")Z = rf" + v" + v'y + m/ 

which is easily integrated to give 



(3.26) 



1 



y' 2 e x 



A / y'e x dx + B + / y'e x (rj" + rj + rjy) dx 



(3.27) 



In addition to the three symmetries of ( |3.24| ) we have the general symmetry 

1 



G 



/ { ^ / y ' eX ( 7? " + r] ' + W) &c\d x + rjd y (3.28) 



for arbitrary functions r\. If we apply the third extension of fl3.28|) to the equation ( p.20[) , 
we obtain, after some simplification, 



Tjy'ti + rj'h = {rjy)' - rf , 



(3.29) 



in which the prime on rj represents differentiation with respect to x and the prime on h(y) 
differentiation with respect to y. If we make use of the chain rule, we may write ( 3.2S| ) as 



( v h)' = {rjy)' - rj' , 



(3.30) 



where now all primes represent differentiation with respect to y and it is evident that rj 
must be a function of y only. We integrate ( |3.30D to obtain 



h{y) = - + y - 1 
V 



(3.31) 



and so, up to a constant, we obtain (2.1) if we take rj to be a constant. 
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4 Reconstruction of the model equation 

The third order equation 

y>» + y» + Ky> + yy> = Q (4.1) 

has the four symmetries 



Gi = 


= «■> G 2 - ^J y/2ex j d x , 


G 3 = 


^(- e lH) 4 ' 


G 4 = 


/■f 1 


/ yy'e x dx 


\d x + d y 


/ W'V 



(4.2) 

which specify it completely. This, of course, is not the model equation with which we 
started and we would do well to verify that the presence of the additional term makes no 
difference to the behaviour of the solution of the equation. 

When we perform the Painleve analysis we find that the first and fourth terms in ( t4.l|) 
are dominant, the singularity is a double pole and the resonances are r = —1,4 and 6 as 
for the model equation, (p.l|). When we make the substitution 



y = Y, a ^~ 2 ( 43 ) 



i=0 

to check for compatibility, we find that 



a = -12, oi = — , a 2 = — - K, a 3 = — , (4.4) 

but then, at the resonance r = 4, a^ = and so there can be no Painleve property for 



(4.1 



Equation (fLW) is trivially integrated to give 



y" + y' + Ky+^y 2 -I = 0, (4.5) 



where / is the value of the first integral. When we perform the Painleve analysis on fl4.5| ), 
we find that the singularity is a double pole and the resonances are at r = — 1 and 6. We 
make the same substitution as in ((O^) and find that the coefficients are given by 

12 1 1 

a = -12, oi = — , a 2 = — - A, a 3 = — — , 

5 25 125 

7 21 + K 2 158 11(2/ + (K + l) 2 ) 

a 4 = , a 5 = — , (4.6) 

2500 20 25000 300 v ; 

that a6 is arbitrary and that there is compatibility at this resonance provided 

2I + K 2 = -—. (4.7) 

625 y ' 
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Finally we compute the Lie point symmetries of fll.q ) when the constants are related 
according to (|Q|). They are 



Gi = & 



and 



Ch 



a */5 



9x + 



6 

25 



K-y)d, 



(4.8) 



We observe that these features are identical to those of ( |2.l[) when we make the substitution 
K = 0. It is not feasible to make numerical experiments for arbitrary K, but we are 
confident that the behaviour of the solutions of fl4.1| ) is the same as the solutions of our 
original model equation. 

We can demonstrate a basis for that confidence by using the symmetries to reduce the 
order of (4J). Naturally we commence with the only point symmetry, G±, which has the 
invariants u = y and v = y' . Equation Q4.5|) becomes 



vv" + v' 2 + v ' + K + u 
and the other symmetries become 







(4.9) 



G 



G 3 ^ 

Ga — > d u 







'\ 


" 


I' du 


}\ 






h -> 1 - exp 
\v 


-J- 


J dv, 


exp 


["/v 


/exp 


/>: 


y 


1 

- ex] 

V 


3 


[-/' 


u 


Jue: 


<;p 


[/t 


du 



a. 



(4.10) 



We observe that G2 has become an exponential nonlocal symmetry and consequently can 



be used to reduce the order of (R 



We obtain the invariants z = u and w = vv' + v 



and ( |4.9p is reduced to the first order equation 

-{K + z). 



dw 
d7 



(4.11) 



We can now see that the presence of K has no bearing upon the behaviour of the solution 
of the equation except as a translation of the initial condition on y(xo), thereby justifying 
our comment above. Under this reduction of order the two unused symmetries become 



G 3 ^d u 



Ga — > zd v 



(4.12) 



and so we can infer that these two symmetries have the Abelian algebra 2A\. 
As (4.11) has two Lie point symmetries, it is obviously integrable. We obtain 

w (z) = C- l -(K + z? 

so that the next equation to be solved is 

vv' + v = C [K + u) 2 . 



(4.13) 



(4.14) 



Equation ( 4.14] ) is an Abel's equation of the second kind which, not surprisingly, does not 
yield to any of the methods of solution proposed by Kamke |11, p. 24, Art 4.10 and 4.11]. 
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In the integrable case the second order equation to be solved is 

i 1 8 

y" + y' + -(y + K) 2 - — = 0. (4.15) 

This can be transformed to an equation of Emden-Fowler type 

w + w 2 = (4.16) 

by means of the Kummer-Liouville transformation 

y = 2e 2x / 5 w(t) - K - — , t'(x) = e"^ 5 . (4.17) 



The solution of ( [4.16| ) is in terms of elliptic functions and so the solution of ( |4.15 ) is in 



terms of elliptic functions of somewhat complicated argument. It is interesting that on the 
surface in the space of initial conditions specified by the value of the integral this solution 
is analytic except for the polelike singularity and yet just off this surface one may find 
nonintegrable behaviour. 

We can always make use of a Kummer-Liouville transformation to bring (|4.5| ) into the 
form of a generalised Emden-Fowler equation. If we put 

y + K = u(x)v(t) + a and t = t(x), (4.18) 

becomes 



{ut' 2 ) v + (2u't' + ut" + ut') v 
+ (u" + u+ua)v + -u 2 v 2 + -a 2 - -j3 2 = 0, (4.19) 

where we have written 2/ + K 2 = f3 2 . We remove the constant term by setting a = j3. We 
remove the coefficient of v by solving the equation 

u" + u' + u(3 = Q. (4.20) 

The coefficient of v is removed if we set t' = u~~ 2 e~ x and we are left with the equation 

v + <p(t)v 2 = 0, (4.21) 

where 

<f>(t(x)) = \e~ x/2 {Ae^ x + Be~^) 5 , (4.22) 



A and B are constants of integration and 7 = J| — (21 + K 2 ) 1 / 2 . We note that the 
form of (p(t(x)) is not one of the functions for which the generalised Emden-Fowler equa- 
tion Q4.21|) , has been found to be integrable [|J, [2JJ . 
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5 Discussion 

As in the case of Krause's analysis for the complete symmetry group of the Kepler problem, 
we can specify our model equation up to a constant. The Kepler problem, being a three- 
dimensional system of second order equations, required eight symmetries, three of which 
were nonlocal. Our scalar third order equation is similarly specified with four symmetries, 
three of which are nonlocal. 



Already in the paper by Krause [12] there was some concern as to the proper way to 
consider the algebraic properties of the symmetries, particularly the nonlocal symmetries. 
In fact he comments "these are not point transformations and one does not obtain from 
them a set of differential operators satisfying a Lie algebra, as in the standard manner." 
]l2| , p. 5738]. The method used by Krause to overcome this obstacle was to look at the 
transformations generated by the symmetries in the coordinates and the rate of change of 
local time. In both of these the nonlocality of the symmetries did not appear and he was 
able to construct a satisfactory group representation. 

Our approach in this matter is somewhat different because we attempt to obtain a 
representation of an algebra from the symmetries under consideration and consequently we 
do not want to have the ambiguities of the algebras of nonlocal symmetries to intrude as an 
unnecessary complicating factor. Our aim is to present the concept in the simplest possible 
context. We also consider the question of the uniqueness of the complete symmetry group, 
because we do not believe that it is unique. Rather one would expect there to be equivalent 
representations and we demonstrate that with the great paradigm, the one-dimensional 
free particle. The free particle in one dimension with the equation 

x = (5.1) 

has the eight Lie point symmetries 

G\ = d x , G2 = t, d x , G3 = xd x , G4 = dt, 

G 5 = 2td t + xd x , G 6 = t 2 d t + txd x , G 7 = xd t , G 8 = txd t + x 2 d x (5.2) 

with the algebra si (3, R) which has the substructure sl(2, R) © s {2A\ © A2 © A{\. If we 
assume the general form of a second order equation, 

x = f(t,x,x), (5.3) 

we find, for example, that / must take the value in ( [5 .1] ) with the combinations of symme- 
tries Gi, G2 and G3; G3, G7 and G&. The algebra of the each triple is ^3,3 (^4> _Dffi s T2, the 
semidirect sum of dilations and translations). The triple G\, G-j and G% also completely 
specifies the equation, but it does not close. To obtain a closed algebra one must add 
the symmetries G3, G4 and G5. The triple G4, G5 and Gq, constituting the subalgebra 
sl(2,R) and characteristic of all differential equations of maximal symmetry, gives the 
Ermakov-Pinney equation 

x = ^ (5-4) 



and needs an additional symmetry, say G\, to recover (p. 
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This example demonstrates quite clearly the possible nonuniqueness of the complete 
symmetry group since these three algebras are quite distinct and cannot be different 
realisations of the same group because they are of different dimension. 

We conclude that we have demonstrated that even an ordinary differential equation 
which has Painleve properties consistent with nonintegrable behaviour is defined in terms 
of its symmetries. We do not believe that our model equation, (|2.lD , is exceptional and 
propose that all ordinary differential equations, be they derived in a Newtonian or rela- 
tivistic context or elsewhere, can be defined in terms of a set of Lie symmetries. We repeat 
the point made earlier at the conclusion of Section 2. The Lie symmetries of a differential 
equation both provide a means to identify the specific structure of the equation and a pos- 
sible route towards its solution. The second property is not automatically available in any 
given symmetry. For a symmetry to be part of the route to the solution of a differential 
equation the symmetry must be able to provide a means for the reduction of order, if not 
immediately, at some level in the reduction process. 
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